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Abstract 

We study the effective energy-momentum tensor (EMT) for cosmological per- 
turbations and formulate the gravitational back-reaction problem in a gauge 
invariant manner. We analyze the explicit expressions for the EMT in the 
cases of scalar metric fluctuations and of gravitational waves and derive the 
resulting equations of state. 

The formalism is applied to investigate the back-reaction effects in chaotic 
inflation. We find that for long wavelength scalar and tensor perturbations, 
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the equation of state is de Sitter-like. 
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I. INTRODUCTION 

It is well known lEHi that gravitational waves propagating in some background space-time 
have an effect on the dynamics of this background. A convenient way to describe the back 
reaction of the fluctuations on the background is in terms of an effective energy-momentum 
tensor (EMT). In the short wave limit, when the typical wavelength of gravitational waves 
is small compared with the curvature of the background space-time, they act as a radiative 
fluid with an equation of state p = \ p (where p and p denote pressure and energy density, 
respectively). The presence of a large amount of gravity waves in the early Universe can 
lead to important consequences for cosmology. For example, it can speed up nucleosynthesis 
and lead to a higher fraction of helium, resulting in constraints on models producing a too 
large amplitude of gravitational waves in the early Universe (see e.g. Ref. Q and references 
quoted therein). 

In most models of the early Universe, scalar-type metric perturbations are more impor- 
tant than gravity waves. On length scales smaller than the Hubble radius, the amplitude 
of scalar fluctuations increases in time, and in most models, scalar perturbations are re- 
sponsible for seeding structure in the Universe. In this paper, we study the back reaction 
problem for both scalar and tensor perturbations (cosmological perturbations and gravita- 
tional waves, respectively). We derive the effective EMT which describes the back- reaction 
and apply the result to calculate EMT for both long and short wavelength fluctuations in 
particular models for the evolution of the Universe. 

One of the main puzzles to be solved is the problem of gauge invariance of the effective 
EMT. As is well known (see e.g. Ref. for a comprehensive review), cosmological pertur- 
bations transform non-trivially under coordinate transformations (gauge transformations). 
However, the answer to the question "how important are perturbations for the evolution 
of a background" must be independent of the choice of gauge, and hence the back-reaction 
problem must be formulated in a gauge invariant way. 

In a recent Letter we demonstrated how the back reaction problem can be set up in a 
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gauge invariant manner. We applied the result to estimate the magnitude of back reaction 
effects in the chaotic inflationary Universe scenario. In this paper, we study in more detail 
the effective EMT of cosmological perturbations. In particular, we derive the equation of 
state satisfied by this EMT. As we show, the back reaction effects of gravity waves and of 
scalar fluctuations decouple. In the short wavelength limit, we recover the result p = \ p for 
gravity waves. In the long wavelength limit, scalar fluctuations about a de Sitter background 
have an equation of state p ~ —p with p < 0. 

The study of back-reaction effects for gravitational waves goes back a long way. Following 
pioneering work of Brill and Hartle , Isaacson 'i defined an effective EMT for gravitational 
waves which was shown to be gauge-invariant for high-frequency waves after averaging over 
both space and time. This prescription only makes sense, however, when considering fluc- 
tuations on scales much smaller than those characterizing the background. In applications 
to physics of the very early Universe the fluctuations of interest have wavelengths larger 
than the Hubble radius and a frequency smaller than the expansion rate. Hence, Isaacson's 
procedure for defining an effective EMT is inapplicable. 

Back reaction effects for density inhomogeneities have been considered only recently, and 
even then without addressing questions of gauge-dependence. The focus of the early work 
of Futamase et al @ and of Seljak and Hui was on effects of inhomogeneities on local 
observables such as the expansion rate of the Universe. For a recent study of this issue 
in the context of Newtonian cosmology, the reader is referred to the work of Buchert and 
Ehlers The focus of our work, on the other hand, is to formulate the back-reaction 
problem in General Relativity in a gauge-invariant manner. 

The outline of this paper is as follows. In the following section we formulate some useful 
properties of diffeomorphism transformations. The back reaction problem is set up in Section 
3 and then in the next section we recast the back reaction problem in terms of gauge invariant 
variables. In Section 5 we first demonstrate that the contributions of scalar and tensor 
fluctuations to the effective EMT do not interfere in the leading approximation. We then 
study in detail the effective EMT for scalar perturbations, focusing on the equations of state 

3 



obtained in the long and short wavelength limits. In Section 6 we derive the effective EMT 
for gravitational waves. As an application, we consider the back reaction of cosmological 
perturbations in the chaotic inflationary Universe scenario. We summarize our results in 
Section 8. 



II. GAUGE TRANSFORMATIONS 

The gauge group of general relativity is the group of diffeomorphisms. A diffeomorphism 
corresponds to a differentiable coordinate transformation. The coordinate transformation 
on the manifold Ai can be considered as generated by a smooth vector field C,"". Let us 
take some coordinate system on Ai in which for instance some arbitrary point P of that 
manifold has coordinates Xp. The solution of the differential equation 



dX 



r [x{P; A)] , (1) 



with initial conditions 



;^"(P;A = 0) = X^ (2) 

defines the parametrized integral curve x"(A) = x'^iP] A) with the tangent vector ^"(xp) 
at P. Therefore, given the vector field on Ai we can define an associated coordinate 
transformation on At as, for instance, Xp-^Xp = A = 1) for any given P. Assuming 

that is small one can use the perturbative expansion for the solution of equation (|l|) to 
obtain 

x'^P = A = 1) = + Cixp) + + 0{e) , (3) 

which we can write in short-hand notation as 

X^(P;A = l) = (e«'if?x")p (4) 

Thus the general coordinate transformation x -^^ x on generated by the vector field ^"(x) 
can be written as 



a 



^ f " = i^x" (5) 

Conversely, given any two coordinate systems and on M. which are not too distant, 
we can find the vector field ^"(x) which generates the coordinate transformations x — *■ 
X in the sense (|^). Of course, in such a way we can not cover all possible coordinate 
transformations 0. However, the class of transformations described above is wide enough 
for our purposes. 

The variables which describe physics on the manifold M. are tensor fields Q. Under 
coordinate transformations x ^ x the value of Q at the given point P of the manifold 
transforms according to the well known law 

Note that both sides in this expression refer to the same point of manifold which has different 
coordinate values in different coordinate frames, that is xp 7^ xp. The question about the 
transformation law for the tensor field Q can be formulated in a different way. Namely, given 
two different points V and P, which have the same coordinate values in different coordinate 
frames, that is xp = x-p, we could ask how to express the components of the tensor in the 
coordinate frame x at the point V (denoted by Qf,) in terms of Q and its derivatives given 
in the frame x at the point V. The answer to this question is found with the help of Lie 
derivatives with respect to the vector field ^ which generates the appropriate coordinate 
transformation x x according to (|^). Its infinitesimal form is given in some books on 
General Relativity (see, e.g. Ref. fl^ ): 

Q{xp = xo) = [Q- C^Q + OieMxv = xo) , (7) 

where denotes the Lie derivative with respect to the vector field C,- Transformation (|^) is 
the infinitesimal form of the gauge transformations of the diffeomorphism group. The finite 
form of it is obtained by exponentiating (^: 

Qix) ^ Qix) = (e-^«g)(a;) 

= Q{x) - {C^Q){x) + ^-{C^C^Q){x) + 0{e) (8) 



Equations (J^) and (P) are tensor equations, where for notational convenience, tensor indices 
have been omitted. Note that, despite of the fact that the transformation law is the 
consequence of transformation law (^, they are different in the following respect: transfor- 
mation law d^) is well defined for any tensor given only at the point P, while (|]) is defined 
only for tensor fields. 

In the derivations which follow below we make substantial use of some properties of the 
Lie derivatives and elements of functional calculus. Therefore for the convenience of the 
reader we would like to recall some basic useful facts from functional analysis and from the 
theory of Lie derivatives. Readers not interested in these formal considerations may skip to 
the next section. 

Let us consider a tensor field G{x) (e.g., Riemann tensor or Einstein tensor) which is 
formed from the metric tensor g{x) and its derivatives dg/dx, that is 

G{x) = G[d/dx,g{x)] . (9) 

Applying the operator exp(— £g) to G{x), we obtain the value of this tensor, denoted G{x) at 
the point V, whose coordinates in the new coordinate frame are x{x) = x. On the other hand, 
G{x) in the frame x can also be calculated from the metric tensor g{x) = {exp{—C^)g){x) and 
its derivatives dg{x)/dx, according to the prescription @. The results should coincide 
and therefore one can conclude (see Ref. [0 for an explicit proof of this nontrivial fact) 
that 



{e-^^G){x) =G 



^,{e-^^g)ix) 



(10) 



that is the Lie derivative can be taken though the derivatives d/dx without "changing" them 
in expressions where these derivatives are used to build the tensors (e.g., Riemann tensor) 
out of the other tensors (e.g., metric tensor). 

We will be interpreting functions G{x) = G[d/dx, g{x)] defined on the manifold Ai as 
the parametrized set of functionals Gx defined on the space of functions g{x') according to 
the formula 
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= G[d/dx, g{x)] = J G[d/dx', g{x')]6{x - x')dx' , (11) 

where 6{x — x') is the Dirac delta function. Then the functional derivative 6Gx/Sg{x') can 
be defined in the standard way 

SG, = I {5GJ5g{x'))5g{x')dx\ (12) 

where 5Gx is the change of the functional Gx under an infinitesimal variation of g{x'): 
g{x') g{x') + 5g{x'). 

If, for instance, Gx = g{x), then 

SGx = Sg{x) = j 5{x- x')Sg{x')dx' , (13) 

and comparing this formula with (|12D we deduce that 

= S(x - x') . (14) 



6g{x' 

As another example, consider Gx = d'^g{x)/dx'^. Using the definitions ([TI|) and (0) one 
gets 

In the following, the functional derivative Fxx' = ^Gx/Sg{x') will be treated as an oper- 
ator which acts on the function f{x') according to the rule 

Fxx'* fix') = J Fxx'f{x')dx' (16) 



We will also use DeWitt's condensed notation |]T6[ and assume that continuous variables 
(t, X*) are included in the indexes, e.g., A°(x*,t) = A*^"'^''*^ = A", where a is used as the 
collective variable to denote {a,x\t). In addition we adopt as a natural extension of the 
Einstein summation rule that "summation" over repeated indexes also includes integration 
over appropriate continuous variables, e.g., 

A'^Ba = f A-{x, t)B^{x, t)dxdt (17) 



We shall write functional derivatives using the following short hand notation: 



6g{x') 6g 

For instance, the useful formula 



/:,G{x) = J d'x'^^^C,g{x') (19) 
which follows from (|lO|) , in condensed notation takes the form 

C^G = G,a{C^gY , (20) 
where in addition we omitted all "irrelevant" indexes. 

III. BACK REACTION PROBLEM FOR COSMOLOGICAL PERTURBATIONS 

We consider a homogeneous, isotropic Universe with small perturbations. This means we 
can find a coordinate system [t, x*) in which the metric (gfiu) and matter (</?) fields, denoted 
for brevity by the collective variable g", can be written as 

q'' = qo + 5q\ (21) 

where the background field Qq is defined as homogeneous part of on the hypersurfaces 
of constant time t and therefore depends only on the time variable t (we recall that the 
variables t and are included in the index a) . The perturbations depend on both time 
and spatial coordinates, and by assumption they are small: 

m « go" • (22) 

From our definition of the background component it follows that the spatial average of 
6q°' vanishes: 

(5g«)= hm ^^^^ = 0. (23) 
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Spatial averaging is defined with respect to the background metric, not with respect to the 
perturbed metric as was done in Ref. [|10|. Our definition is the appropriate one when 
estabhshing what "perturbations" are and when constructing the general back reaction 



framework. The averaging of Ref. |]T0[ is appropriate when discussing the "expected" values 
of physical quantities for real observers in systems with fiuct nations on scales smaller than 
the Hubble radius. 

The Einstein equations 

G^, - StiGT^, := = (24) 

can be expanded in a functional power series in Sq"" about the background ^q, if we treat 
G^i, and as functional of q'^ 

n(go") + n,,|^. + ^ n,,,|^g Sq'^Sq' + 0{Sq') = . (25) 

To lowest order the background q^ should obey the Einstein equations 

n(go") = (26) 

and the fiuctuations satisfy the hnearized Einstein equations 

n,a(go)5g'' = o . (27) 



By definition, the spatial average of the linear term in Sq"" in (|25|) vanishes. The term 
quadratic in however, does not. Therefore the spatial averaging of (^) leads to higher 
order corrections in the equations describing the behaviour of the homogeneous background 
mode. Thus the "corrected" equations which take into account the back reaction of small 
perturbations on the evolution of the background are 

n(go")= -i(n,,,5g"5g^) (28) 



At first sight, it seems natural to identify the quantity on the right hand side of 
as the effective EMT of perturbations which describes the back reaction of fiuctuations 
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on the homogeneous background. However, it is not a gauge invariant expression and, for 
instance, does not vanish for "metric perturbations" induced by a coordinate transformation 
in Minkowski space-time. 

In the next section we will rewrite the back-reaction equations in a manifestly gauge- 
invariant form. First, however, we want to show that the physical content of ( pS] ) is inde- 
pendent of the gauge chosen to do the calculation in, provided that we take into account 
that the background variables change to second order under a gauge transformation. 

The coordinate transformation induces (to second order in perturbation variables) 
the following diffeomorphism transformation of a variable q 

q = qo + 6q ^ e~^« (go + Sq) 
= qQ + 6q- C^qo - C^6q+^ Cjqo (29) 

Hence, to linear order, the change in 6q is 

Sq^ Sq = Sq- C^q^ , (30) 

while, to second order, the background variable transforms nontrivially as 

go go = go - + ^(^l^o) , (31) 

where (0 = has been assumed. 

In order to prove that the equations ( PB] ) are independent of gauge, we must show that 

H(go) = -^(n ,fe5g'^ 6q') ^ H (go) = -^(n.^.^r^g') (32) 



(to second order in perturbation variables). Making use of (|30| ) and (|3lD we obtain 

1 



and 



H (go) = n(go) - H,,(£^5g'^) + -H,,(£|g«) (33) 



^ (n,„, SrSq' ) = - ^ (H,,, Sq'^Sq") (34) 
+ {UatC^q^o^q") - liUatC^q'o^^q',) . 
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In order to show that the extra terms on the right hand sides of ( p3D and (p^ cancel out in 



(0), we make use of (|T9[) in the equation 



l)n(go + 5g)) = 0. 



(35) 



Expanding this equation to second order in perturbations, we obtain the identity 



(36) 



which completes the proof that the extra terms mentioned above cancel out. 

IV. GAUGE INVARIANT FORM OF THE BACK REACTION EQUATIONS 

Although we have shown in the previous section that the physical content of Eqs. (^Sf ) 
should be the same in all coordinate systems, it is useful for many purposes to recast these 
equations in an explicitly gauge-invariant way. In particular, this will allow us to define a 
gauge invariant EMT for cosmological perturbations. 

We start by writing down the metric for a perturbed spatially flat FRW Universe 



where a{t) is the scale factor, and where the 3-scalars (j), ip, B and E characterize scalar metric 
perturbations. The symbols Si and Fi are transverse 3- vectors (giving vector fluctuations), 
and hij is a traceless transverse 3-tensor (gravity waves). 

We will now attempt to construct gauge-invariant variables Q from the gauge-dependent 
quantities q. To do that let us take the set of 4 quantities (not necessarily a four-vector) 
and form a Lie operator with (denoted by Cx)-, treating formally as a four- vector. 
Later on, after we specify the properties which should satisfy if we want it to help build 
gauge invariant variables, we will then construct X'^ explicitly out of the metric perturbation 
variables. However, let us leave X^ unspecified for a moment. Using Cx we define the new 
variable Q according to the prescription 



ds 



,2 



(1 + 2(t))dt^ - 2a{t){B^i - Si)dx'dt 



(37) 



a^{t)[{l - 2i,)5ij + 2Ey + + Fj,, + h,j\dx'dx^ , 
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g = e^^g . (38) 

If we attempt the same construction after performing a gauge transformation (with param- 
eter 1^), then taking into account that as a result of this transformation X X and 

q^q = e-^^q , (39) 

we obtain 

Q[x) = e^^e-^iq{x) . (40) 



If we demand that Q is gauge- invariant, that is Q = Q, then comparing ( pS]) and (|40|) we 
arrive at the condition 

= e^^e^« . (41) 
which imposes strong restrictions on the transformation law X X. From (ffTF) and by using 



the Baker-Campbell-Hausdorff formula for the products of exponentials of operators, one can 
easily find that the condition for gauge-invariance of Q implies that under a diffeomorphism 
generated by ^, 

X^-^X^' = X^' + C + -[X,^f + ■■■, (42) 

where ■ ■ ■ denote terms of cubic and higher order. 

Note that Q is a. gauge invariant variable characterizing both the background 

Q-, = q-, + {Cx6q'') + l{Clq^,) (43) 

and the linearized perturbations 

= 6q^ + Cxql (44) 

Making use of the gauge invariant variable Q we can recast the back-reaction problem 
(P5|) in a manifestly gauge invariant form. If q satisfies the Einstein equations, it follows 
from our basic identity (0) that 
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e^^n(g) = n(e^^g) = U{Q) = 0. (45) 

Expanding the above equation to second order in SQ and taking the spatial average of the 
result yields 

U{Qo) = -^{UabSQ^SQ") = —{G^abSQ^SQ') + ATiGiTabSQ^SQ'), (46) 

which is the desired gauge invariant form of the back-reaction equation. Reinserting tensor 
indices, the above equation can be rewritten as 

GM) = SnGlT^^iQo) + r^.(5Q)], (47) 

where 

r,uiSQ) = --J^{n,,,ab5Q''5Q') (48) 

iOTTG 

can be interpreted as the gauge-invariant effective EMT for cosmological perturbations. 
At this point, however, we must return to the question of what is. It is a question 



of linear algebra to find the linear combinations of the perturbation variables of (|37D that 
have (to linear order) the required transformation properties, namely 

Xi^ ^ Xf^ + ^f". (49) 

The solution we will use is 

X^' = [a{B - aE) , -E^ - F^] , (50) 

where a dot denotes a derivative with respect to the time variable t. But this choice is not 
unique. There is a four parameter family of possible choices labelled by real parameters 
a, a, 7 and a. The component is 

X° = aa(5-aE) + (l-a)(-^), (51) 

a 

the X* components have a traceless piece 
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X^, = -aF, + (1 - a) f dt'-Si{t') 

Jo a 



(52) 



and a trace 



XI = X. 



L,i ) 



(53) 



with the function Xl defined as 

Xl = -jE + (1 - 7) /* dt'- \-X\a; t') - B{t'] 

Jo a la 



(54) 



Finally, 



X' = X'^r + X 



(55) 



Demanding regularity in the limit where the expansion rate vanishes forces a = a = 1. In 
this case, the dependence on 7 drops out of Eq. (|5^) and we are left with a one-parameter 
degeneracy of labelled by a. 



V. ENERGY-MOMENTUM TENSOR FOR SCALAR PERTURBATIONS 

In this and the following section, we calculate the effective EMT for scalar and tensor 
perturbations, respectively. Since vector modes decay in an expanding Universe, we shall in 
the following take them to be absent. 

In models such as the inflationary Universe scenario, scalar and tensor modes are statis- 
tically independent Gaussian random fields. In this case, the effective EMT (4.10) separates 
into two independent pieces, the first due to the scalar perturbations, the second due to the 
tensor modes. 

r,um = r^f'^'^m + <™(5Q) (56) 

Note that if we neglect vector perturbations then it is enough to consider only scalar 
modes in X^^. Hence, the contribution to the effective EMT (^S]) coming from the term 
CxQq appearing in SQ"" is a contribution to t^^"''^'" alone. 
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It is easy to verify that for our choice of X'^, namely, 

Xf" = [a{B -aE),-E,i], (57) 

the variables SQ"^ of (0) for the metric perturbations correspond to Bardeen's gauge invari- 
ant variables 0. In fact, the application of (0) yields the following gauge invariant metric 
tensor: 

Sglf = S9,u + Cxgf} , (58) 

and from the time-time and diagonal spatial components, respectively, we immediately ob- 
tain the two gauge invariant variables, 

0(30 ^ $ = + _ Q^)]-^ 

i)^^'^ = m = i) -a{B -aE). (59) 

For the remaining components of the metric tensor, the gauge invariant combinations vanish. 
Hence, calculating the general gauge-invariant effective EMT (^Sp reduces to calculating 

1 



<n^,,„,5g'^5g^> (60) 



in longitudinal gauge (S = £" = 0), in which 

ds^ = (1 + 20) dt^ - a^it) (1 - 2%l)) 6ijdx' dx^ (61) 

For many types of matter (scalar fields included), Tij is diagonal to linear order in Sq. 
In this case, it follows from the linearized Einstein equations that 

(f) = 'ip (62) 

Thus, in longitudinal gauge the variable entirely characterizes the metric perturbations. 
We shall consider scalar field matter, in which case the linear matter fluctuations are de- 
scribed by 6ip. The motivation for our choice of matter follows since we have applications 
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of our formalism to inflationary cosmology in mind. The linearized Einstein equations also 
relate 5(p and (j). In fact, there is a single gauge invariant variable characterizing linearized 
scalar fluctuations. 

We will now calculate t'^J for scalar perturbations (metric ([6T|)). The contribution of 
gravitational part to the effective EMT can be easily calculated with the help of the formulae 
(see, e.g., in pg.965, Eq. (35. 5g) a& b): 



[^9MP ~ ^9i3t,\a) - (63) 



if we substitute in these expressions the metric (^1]). In this formula | denotes covariant 
derivatives with respect to the background metric. Expanding the energy-momentum tensor 
for a scalar field 



(64) 



to second order in 5ip and 5g and combining it with the result for G^^^ we obtain from (^Sp 
roo = ^ [+l2Hm - 3((0)2) + 9a-2((V0)2) 

+ i\/"(v9o)(V) + 2V'{^o){(pS^) , (65) 

and 



{24H'^ + 16H){(j)'^) + 24H 
4 



+ ((0)2)+ 4(00) --a-2((V0)^) 



+ 4vio'(r) 



- -V"{^o){6^') + 2V'M{<p6y^) 



where H is the Hubble expansion rate, and where we have used the fact that 
theories in which 6Tij is diagonal at linear order. 



(66) 

ip for 
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Before discussing the long and short wavelength limits of the equation of state satisfied 
by the cosmological perturbations, let us briefly recall a few crucial points from the theory 
of linear fluctuations (see Ref. for a comprehensive overview). 

The simplest way to derive the equations of motion satisfied by the gauge invariant 
perturbation variables is to go to longitudinal gauge {B = E = 0), in which the gauge 
invariant variables $ and \I' coincide with the metric fluctuations and ip. The equations 
of motion (^) 



derived in this gauge coincide with the equations for Bardeen's gauge invariant variables. 

For scalar field matter, as mentioned above, (p = ip. In this case, the 00 and ii perturba- 
tion equations combine into a second order differential equation for $ which on scales larger 
than the Hubble radius and for a time-independent background equation of state has the 
solution 



(modulo the decaying mode). The Oi equations give a constraint relating (p and 6ip: 



Our aim is to work out r^i, in an inflationary Universe. In most models of inflation, 
exponential expansion of the Universe results because (fo is rolling slowly, i.e. 



5G^ = 8nG5Tlp 



(67) 




(68) 



(p + H(p = AirGipo 6ip . 



(69) 



V 



(70) 



3H ' 



where a prime denotes the derivative with respect to the scalar matter field. The (p term 



in Eq.(p9|) for the nondecaying mode of perturbations is proportional to a small slow roll 
parameter and therefore can be neglected. Thus from Eqs. (|69D and (^) we obtain 



(71) 
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When considering the contributions of long wavelength fluctuations to t^u, we can neglect 
all terms in ( |65D and ( p6| ) containing gradients and factors. Because of the "slow rolling" 
condition ([7^), the terms proportional to 0q and H are negligible during inflation (but they 
become important at the end of inflation). Hence, in this approximation 

Too -^V" < Sip'^ > +2V' < ^Sifi > (72) 

and 

Tij ~ a^6^j l^H^ < > -Iv" < > +2V' <4>S^>} . (73) 
I ttG 2 J 

Making use of (|7l|), this yields 

P^^^o = (^2-^-4vj <02> (74) 

and 

Ps = -^r; = -p(^) . (75) 

Thus, we have shown that the long wavelength perturbations in an inflationary Universe 
have the same equation of state = —p^ as the background. 

One of the main results which emerges from our analysis is that 

Ps < (76) 

for the long wavelength cosmological perturbations in all realistic inflationary models. Thus, 
the effective t^^ counteracts the cosmological constant driving inflation. Note that the same 
sign of p emerges when considering the vacuum state EMT of a scalar fleld in a flxed 
background de Sitter space-time (see e.g. Refs. [jl8|,|19|). 



For short wavelength fluctuations [k ^ aH), both cj) and 5(p oscillate with a frequency 
oc /c. In this case (see (|69D): 

id 

~ At:G — ipoSip . (77) 



Hence, it follows by inspection that all terms containing in (|65|) and (|66|) are suppressed 
by powers of Ha/k compared to the terms without dependence on (p, and therefore 

Too^l< > +\a-^ < {V6^f > +^ V'i^o) < V > (78) 

and 

r^J = a% [I < {6^f > -^a-' < {V6^Y > ~ V"{^o) < V >} (79) 
which is a familiar result. 



VI. ENERGY-MOMENTUM TENSOR FOR GRAVITATIONAL WAVES 



In the case of gravitational waves the metric takes the form 

ds^ = dt^ - a^{t){5ik + hk)dx'dx^ (80) 

where hi^ is defined as the transverse traceless part of the metric perturbations, and therefore 
the components 5gik = hij are gauge invariant themselves. 

In the absence of matter fluctuations, the effective EMT t^^ is given by the first term in 
(^). Making use of the second variation of R^^ given by (0) and the equation of motion 
for gravitational waves 

hij +3^h^, -^V^hij = 0, (81) 



we obtain 



d ■ 1 / • ■ 1 , 

iirGroQ = - <hke hke > +- <hkehke> H — 5- < hke,mhke,m > (82) 



and 



f 3 3 ■ • 1 

SnGTij = 6ija < hke,mhke,ni > —- <hkihki>\ 

+ < hikhkj> +^ < hke,ihke,j > —- < hik/hjk,i > ■ (83) 
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The first expression can be interpreted as the effective energy density of gravitational waves 

Pa^ = ■ (84) 

The relation between Tij and the quantity which we could naturally interpret as an 
effective pressure is not so straightforward as it looks at first glance. The problem is that the 
energy momentum tensor for the gravity waves r^^, is not conserved itself, that is, r^^ 7^ 0. 
This is not surprising since the gravitational perturbations "interact" with the background 
and only the total EMT must be conserved, 

(TOo) + r,^W))|M = (85) 

as a consequence of Bianchi identities for background, G'j^{Qo)\^ = 0. In fact expanding the 
exact conservation law T^^^ = in perturbations to second order and averaging the resulting 
equation we obtain 

T^iQoh = -(^'^rg^) T2 + e^n,) (86) 

Deriving Eq. (^) we took into account that in our case (when we have only gravity waves) 
matter perturbations are absent. The energy momentum tensor for the background is diag- 
onal and isotropic, that is, Tg = p^^^ and = —p^^^Sl- Also, we will consider only isotropic 
fields of gravitational waves. In such a case the conservation law (^), taking into account 
(^), can be written down explicitly in a familiar form as 

Pgnj + SH{Pgy, + P g^) = (87) 

where 

P.. = -^r;-^f)0(p^°^+P^°^) (88) 

can be interpreted as the pressure of gravitational waves 

The second term in (^) does not contribute to the pressure only in a de Sitter Universe 
in which p^^^ = —p^^\ In this case the term —^t- itself can be interpreted as a pressure. 
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A similar but a bit more complicated analysis can be done for the scalar perturbations, for 
which one third of the trace of the spatial part of the EMT can be interpreted as a pressure 
also only in a de Sitter Universe. 

As for scalar perturbations, we will now study the equation of state for gravity waves 
both in the short and long wavelength limits in various models for the evolution of the 
Universe. First, assuming that the field of gravity waves is isotropic and by averaging the 
diagonal elements of (^) we obtain the contribution of r to the pressure 

3~ ^* ~ 24q2 ^ '^ke,miike,m > — <hktiikt> ■ (89) 

For fluctuations with wavelength smaller than the Hubble radius, the first term on the right 

hand side of ( P^ ) and the second term in ( |55D are negligible and the time average of the 

temporal and spatial gradient terms are the same. Hence 

111 

Pgw — gPgui — ^^^2 ^ hke^mhk£,m ^ (90) 

where <C ^ indicates that in addition to spatial average, a time average over a period 
T ^ H^^ has been taken. As expected, short wavelength gravitational waves behave like 
radiation, independent of the evolution of the background, and their energy density decays 
as a~'^(t). 

In the case of long wavelength gravitational waves the calculations are less straightfor- 
ward. First, we consider a de Sitter background : 

a{t) = e^(*-*°) (91) 

where to is a reference time. Let us take an isotropic field of gravity waves with comoving 
wavenumbers k. For the nondecaying mode of long wavelength gravity waves {k/a ^ H), 



the solution of Equation (|8lD is 



i + ^f-iV + o((-if) 



e"""^ , (92) 



2 \aH J ''aH' 

where Ak is a constant (related to the spectrum of gravity waves) and eij is the polarization 
tensor. Substituting this solution in Formula (^) we obtain the following expression for the 
energy density of the gravity waves to lowest order in k : 
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and correspondingly from (^) we derive that the pressure is: 

1 



Pgw 



' Pgw 



(94) 



Note that the whole contribution to the pressure in this case comes from the first term in 
(^) . The energy density of long wavelength gravity waves decays as ~ a~^. 
In a radiation dominated Universe the scale factor increases as 



.to. 



(95) 



and the solution for long wavelength gravity waves {k <^ Ha) is 



+o((— )3) 

6\aH ^aH^ ' 



ikx 



Inserting these results into (^) and (|83D yields 

Pgw - 



1 ^ 



SttG 24 a2 



(96) 



and 



Pgw 



' Pgw 



(97) 



In this case both of the terms in (^8|) give comparable contributions to the pressure, namely 



the contribution of the first term there is pi = —Pgw and the contribution of the second term 



is negative and equal to p2 = —^Pgw As in the previous case the energy density decays as 

^-2 



VII. BACK-REACTION IN INFLATIONARY UNIVERSE 

As an application of the formalism developed in this paper, we will study the effects of 
back-reaction in inflationary cosmology. To be speciflc, we consider a single fleld chaotic 
inflation model @ and take the inflaton potential to be 
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V(^) = lmV- (98) 

Furthermore, we specify an initial state at a time ti in which the homogeneous inflaton field 
has the value ^o{ti) and the fiuctuations are minimal. 

We will focus on the contribution of long wavelength modes to the effective EMT, which 
by assumption vanishes at the initial time ti. Provided that (po{ti) is sufficiently large, the 
slow-rolling conditions are satisfied and exponential expansion will commence. At this point, 
quantum vacuum fiuctuations begin to generate perturbations on scales which "leave" the 
Hubble radius (i.e. whose physical wavelength becomes larger than the Hubble radius). 
As time proceeds, more modes leave the Hubble radius, and hence the contribution to the 
effective EMT builds up. We wish to estimate the magnitude of the resulting effective EMT. 

As discussed in Setion 5, scalar metric perturbations in this model are characterized by 
a single function (pk (where k stands for the comoving wave number). From the constraint 



equation ([7l|) it follows that 

^<^k - -<^o<Pk ■ (99) 



Hence, the dominant terms in the effective energy- momentum tensor T^,y (see (74) and (|75|)) 
are proportional to the correlator (0^). The amplitudes of (pk are known from the theory 
of linear cosmological perturbations. Using the results for 0^ valid during inflation we 
obtain 



dk 2 _ m'^M'^ 

J I rk I e-)r 

k, k 6^ 



f'^t dk 


r H{t)a{t)-\ 


Iki k 


k 



(100) 



where t denotes physical time, and Mp is the Planck mass. The integral runs over all modes 
with scales larger than the Hubble radius, i.e. 

k < kf{t) = H{t)a{t) , (101) 

but smaller than the Hubble radius at the initial time ti, i.e. 

k > ki = H{ti)a{ti) . (102) 
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The infrared cutoff ki is a consequence of our clioice of initial state. 
For the potential V{ip) considered here, the scale factor is given by 

a{t) = a{t^expi^[ipliU) - ifUt)]) . (103) 

The intregral over k in ( |100| ) can be calculated explicitly, giving 



Ml 



3 



(104) 



327rVoW 

Making use of (|7^ , we finally obtain the fractional contribution of scalar perturbations ps 
to the total energy density 



Ps{t) 3 mVo(^ 



(105) 



po 47r M| 

where po is the background energy density of the homogeneous scalar field ipo. In situations 
in which the ratio ( |105| ) becomes of the order 1, back- reaction becomes very important. 

Several consequences can be derived from (|105|) . First of all, back-reaction may lead to 
a shortening of the period of inflation. Without back-reaction, inflation would end when 
foit) ~ Mp (see e.g. |^). Inserting this value into ( p.05|) , one can expect that if 



then back-reaction will become important before the end of inflation and may shorten the 
period of inflation. It is interesting to compare this value with the scale 

which emerges in the scenario of stochastic chaotic inflation ElB as the "self-reproduction" 
scale beyond which quantum fluctuations dominate the evolution of <fo{t). Notice that 
(Psr ^ Vbr since m <C Mp. Hence, even in the case when self-reproduction does not take 
place, back-reaction effects can be very important. 



Alternatively, we can fix (po{ti) and use the expression ( |105| ) to determine at which value 
of (fo (denoted by (fo{tf)) back-reaction becomes important, which presumably implies that 
inflation will end at that point. The result is 
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We will conclude this section with an analysis of the back-reaction of scalar perturbations 
on the evolution of the homogeneous component of scalar field ipo{t)- The equation for the 
scalar field ifo{t) taking into account the backreaction of perturbations can be obtained if 
we start with the exact Klein-Gordon equation 

□30+53 ((^o + M + ^'(<^o + M = , (109) 

expand it to second order in perturbations 5g, Sip and take the average. The result is 

(00 +3if¥3o)(l + 4(02)) + r + iv'"(V) - 2(050) 

-4(0(50) - 6if(0(50) + 40o(00) - ^(0V'(5yp) = . (110) 

For long wavelength perturbations in inflationary Universe, the terms containing spatial or 
time derivatives of and 6ip are negligible. Hence, for the potential (0), the equation (|110|) 
becomes 

2 

d 771 

^0 + 3-00 = -Y^l(^fo . (Ill) 

We conclude that back-reaction of perturbations on the evolution of ipo become very impor- 
tant when (02) ~ 1, at the same time as they become important for the evolution of the 
background geometry of space-time. 



VIII. CONCLUSIONS 



We have defined a gauge invariant effective energy-momentum tensor (EMT) of cosmo- 
logical perturbations which allows us to describe the back-reaction of the perturbations on 
the evolution of background space-time . Our formalism can be applied both to scalar and 
tensor perturbations and applies independent of the wavelength of the fluctuations. 

In particular, our analysis applies to cosmological perturbations produced during an 
period of inflation in the very early Universe. In this case, we have worked out the specific 
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form of the effective EMT for both density perturbations and gravitational waves. The 
contribution of long wavelength scalar fluctuations to the energy density has a negative sign 
and thus counteracts the cosmological constant (note that this effect is a purely classical 
one, in contrast to the quantum mechanical effects counteracting the cosmological constant 
discussed in Refs. |]TB|,0,pB[ ) . The equation of state of the effective EMT is the same as that 



of the background. The contribution of long wavelength gravitational waves to the energy 
density also has a negative sign, and in this case the equation of state is p = — l/3p. Note 
that an instabihty of de Sitter space to long wavelength fluctuations was also discovered in 
Ref. [|4|. 



Applied to the chaotic inflationary Universe scenario, we found that the back-reaction 
of generated perturbations on the evolution of the background can become very important 
before the end of the inflation even if we start at an energy scales below the "self-reproduction 
scale" . 
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